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ABSTRACT 

Recent reports suggest that elliptical galaxies have increased their size dramatically over 
the last ~ 8 Gyr. This result points to a major re-think of the processes dominating the late- 
time evolution of galaxies. In this paper we present the first estimates for the scale sizes of 
brightest cluster galaxies (BCGs) in the redshift range 0.8 < z < 1 .3 from an analysis of deep 
Hubble Space Telescope imaging, comparing to a well matched local sample taken from the 
Local Cluster Substructure Survey at z ~ 0.2. For a small sample of 5 high redshift BCGs 
we measure half-light radii ranging from 14 — 53kpc using de Vaucuoleurs profile fits, with 
an average determined from stacking of 32.1 ± 2.5 kpc compared to a value 43.2 ± 1.0 kpc 
for the low redshift comparison sample. This implies that the scale sizes of BCGs at z = 1 
are ~ 30% smaller than at z = 0.25. Analyses comparing either Sersic or Petrosian radii 
also indicate little or no evolution between the two samples. The detection of only modest 
evolution at most out to z — 1 argues against BCGs having undergone the large increase in 
size reported for massive galaxies since z — 2 and in fact the scale-size evolution of BCGs 
appears closer to that reported for radio galaxies over a similar epoch. We conclude that this 
lack of size evolution, particularly when coupled with recent results on the lack of BCG stellar 
mass evolution, demonstrates that major merging is not an important process in the late time 
evolution of these systems. The homogeneity and maturity of BCGs at z = 1 continues to 
challenge galaxy evolution models. 
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1 INTRODUCTION 

There is currently great interest in the evolution of the scale size of 
massive galaxies over the age of the Universe, as this provides con- 
straints on the merger history and feedback processes that govern 
galaxy evolution. However, this is still a controversial field as there 
have been claims that massive galaxies at high redshif t are smaller 
than galaxies with a s imilar mass i n the local Univer s e {Daddi et aT 
200. 



mained unchanged and still exist today I Valentinuzz i et alJl20TCt 
ISaracco et al1l2010h . 

In this paper we look for evidence of a size evolution in 
brightest cluster galax ies (BCGs), the most luminous and homoge- 



nous class of galaxy JSandage|ll972l 1 19761; ICollins & Mannll 19981 : 



Whilev et al. 2008k Istott et~ 



2008 



Collins et al. 



2009). 



l galaxies with a similar mass in the local Universe (Daddi et al 
'5b iTruiillo etafl 120061, [20071; iToft et all |2007|; ICimatti et al 



2008MFranx et al.l2008l:lvan der Wei et alj2008b|van Dokkum et al 
2008b lDamianov et al.ll2009blvan Dokkum et al.ll2009ll2010l): how- 
ever others reach different conclusions (e.g. iMancini et al .11 20 101) 
and one major criticism is that any systematic errors involved in 
the high redshift observations and a nalysis would tend to bias 
a result towards such an evolution dPriverl l201(j|) . In addition, 
there are observations of massive and compact galaxies in the lo- 
cal Universe and therefore, rather than constituting an evolution, 
the massive compact galaxies seen at high redshift may have re- 
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The ad- 
vantage of studying BCGs in particular, is that they are also the 
largest galaxies and can therefore be seen at great distances; but 
most importantly, they reside at the centres of rich clusters which 
are thought to be comparable environments across all epochs. 

There have been a large number of studies of the surface 
bright ness profiles and scale sizes of BCGs in the local Universe 
(e.g. IOemlerlll976b ICarteilll977b ItTpessel et al.|[l987b ISchombertl 
1 19861 1 19871; lOegerle & Hoessell 1 199 lb iGraham et alj Il996l: 
ii etalj |2007| 



Bernardi et al. 2007; I von der Linden et al l 1 20071 : iFasano et al.l 



2010) with measured scale sizes found to range typically from 
^lOkp c to ^lOOkpc. In terms of evolution in the scale size of 
BCGs, iBernardil d2009h find that when comparing samples in an 
internally consistent way there is a 70% increase in the size of 
BCGs in the 3 Gyr between z — 0.25 an d z ~ 0, thought to b e due 
to dry, dissipationless minor merging. lAscaso et also 
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find an increase in size between z ~ 0.4 and z ~ 0, although with 
no measured change in the shape of the profile, which they claim 
cannot therefore be due to merging and is thus caused by feedback 
processes in the centre of the galaxy. 

Measuring the scale size of la rge galaxies and B CGs in par- 
ticular is notoriously difficult (see lLauer et al J [2007b and one of 
the principle concerns of interpreting the scale-size evolution of 
BCGs and elliptical galaxies is the significant disparity in published 
va lues even at low redshift. Fo r example, the scale radii measured 
bv lvon der Linden et al.l J2007T) are signi ficant ly smaller than t hose 
measured by both iBernardi et al.l d2007t) and iLiu et al.l d2008l) for 
essentially the same Sloan Digital Sky Survey (SDSS) data. At 
moderate redshift there is also significant disagreement. As an il- 
lustration one particular BCG (in the cluster MS 0451—03) with a 
scale size m easured using a de V aucouleurs profile from near-IR 
HST data bv lNelson etafl j2002t) is found to be a factor 6 smaller 
than a similar optical analysis of the same galaxy by Bildf ell et al.l 
d2008l): in th i s case possibly due to the small field of view of the 
iNelson et alj d2002h observation and the presence of a foreground 
spiral galaxy dominating the light at larger BCG-centric radii. In 
other cases these differences appear to be due to the subtly differ- 
ent ways that the scale size of BCGs is defined or measured and also 
the variable depth and quality of the data between samples. Some 
authors fit only de Vaucoleur profiles while others fit free Sersic 
curves, which can give significantly differe nt results, and in s ome 
cases double fits (e.g. Sersic + exponential. I Ascaso et al ■l2010h arc 
used which can be difficult to interpret or compare with any other 
investigations. In addition, some stud ies only fit circularly symmet- 
ric pro files (e.g. iBildfell et alj|2008l) while others include elliptic- 
ity (e.g. lLiu et alj200a) which can give di fferent results depending 
on the galaxy's shape (Nelson et al. 2002). There is also the prob- 
lem t hat Sersic index and scale size are coupled (e.g lGraham et al.l 
fl996h and as such when fitting a free Sersic curve one can find 
answers that are degenerate in these two parameters. In terms of 
the data, although the Sersic fit is parametric and therefore depends 
only on the shape of the surface brightness profile, the depth of the 
data or the limiting surface brightness used for the fitting process 
may in some studies be inadequate to see enough of the profile to 
return a good description of the galaxy's lig ht distribution, shown 
to be an important factor bv lLiu et alj a2008l) when explaining the 
discre pancy between their results and those o flvon der Linden et al.l 
(2007). Different treatments of the background subtraction, partic- 
ularly at longer wavelengths, and the masking out of near neigh- 
bours in crowded fields may also lead to a discrepancy between 
results from different studies. In addition, significant atmospheric 
seeing from ground-based measurements, particularly when going 
to moderate redshifts, if not correctly deconvolved with the galaxy 
image, may cause erroneous fits to the inner parts of galaxies af- 
fecting the entire measured profile. 

BCGs are generally found to lie off the Kormendy relation 
(surface brightness versus scale size, iKormendvl 1 19771) . having 
larger radii for a given surfa ce brightness compared to the general 
elliptical galaxy populati on dHoessel et al.lll987l ; ISchombertlll987l ; 
Oegerle & Hoes selll99lh . A similar relation to this, luminosity ver- 
sus sc ale size, is found to be steeper for BC Gs than for regular ellip- 
ticals jBernardi et al.l2007l;lLiu et al.l2008h . although other authors 
using essentially the same data hav e not seen evidence for a steeper 
slope dvon der Linden et alj|2007t) . The offset of BCGs from the 
Kormendy relation of regular ellipticals is most likely due to the 
extended cD halo associated with BCGs in rich clusters, which 
is a low surface brightness component o f the galaxy light profile 
dRichstone|[l976l ; lLopez-Cruz et al.ll 19971) . The origin of this halo 



is thought to be the addition of tidally stripped mater ial, from other 
cluste r members, to the outer regions of the BCG dHoessel et al.1 
1987'). The halo itself however, provides difficulties for the study 
of BCG sizes as this is sometimes found to be a separate compo- 
nent of the surface brightness, potentially requiring a double fit, 
or it can modif y a single fit to b e more like a power law (Sersic 
index n — > oo, ISchombed[l986l) rather than a de Vaucouleur pro- 
file (Sersic index, n = 4), which therefore does not conver ge and 
can continue into th e intra-cluster light (ICL) for 100s kpc i Carted 
1 19771; IOemledfl976l) . Estimates suggest as much as 33 - 40% of 
an entire cluster's optical lumin osity is contained in the BCG + 
ICL of which 80% is in the ICL (Go nzalez etai]|2007h . To account 
for this many investigators fit BCG light profile s to a consistent 
surface brightness level (e.g. lGraham et al.lll996l) as this has been 
shown to signifi cantly affect the recovered value for the galaxy size 
dLiu et al.ll2008l) 

Our search for evolution in the scale-size of BCGs is timely 
considering our recent work which suggests that BCGs as a popu- 
lation undergo little evolution in the ir stellar mass out to z = 1.5 
dCollins et al.ll2009l ; ls"tott et alj|2010l) . Any evolution we see in the 
scale size would suggest a number of possible processes dominat- 
ing the late-time evolution of BCGs: minor merging; interactions 
which add mass to the outer regions of the galaxy; feedback from 
an active galactic nuclei (AGN); or a central starburst that disrupts 
the mass in the centre of the systems giving rise to adiabatic ex- 
pansion dHopkins et al ]|2010l) . In any case from a theoretical point 
of view, brightest cluster galaxies should be ideal candidates with 
which to study size and shape evolution as a result of hierarchi- 
cal assembly, because they are thought to have undergone many 
mergers and interactions over their history as a result of being lo- 
cated so close to the centre of mass of the largest dark matter haloes 
dDe Lucia & Blaizotll2007l) . 

This paper sets out the first comprehensive study of BCG scale 
sizes in a small X-ray selected sample at high redshift (z > 0.8) 
and compares the results to a similar and larger sample in the local 
Universe (z ~ 0.2), exploiting the best space-based imaging data 
available, well matched in both rest-frame filter and depth, in an 
attempt to reduce the problems which beset other studies discussed 
above. 

A Lambda Cold Dark Matter (ACDM) cosmology (Q M =0.3, 
Qvac =0.7, Ho =70 km s _1 Mpc -1 ) is used throughout this work. 



2 THE SAMPLE 

Table Q] details our sample of 5 of the most distant, spectroscop- 
ically confirmed, galaxy clusters. The sample consists of clusters 
discovered by several X-ray surveys. The clusters all have red- 
shifts in the range 0.8 < z < 1.3 and have X-ray luminosities 
of 1 i$ L x < 18 x 10 44 ergs _1 (0.1 - 2.4 keV). All of these clus- 
ters have been imaged with the Hubble Space Telescope Advanced 
Camera for Surveys (HST/ ACS) through the F850LP filter with the 
exception of CLJ1226 + 3332 which is imaged with the F814W 
filter. We choose these particular 5 clusters as they are the highest 
redshift clusters with the longest image exposure times (> lOks) 
which is key for observations of low surface brightness features that 
may otherwise be difficult to detect due to the strong dependancy 
of cosmological surface brightness dimming on redshift. The data 
for CLJ0152.7- 1357, MS1054.4 - 0321, RDCS J0910 + 5422 
and RDCS J1252.9 - 2927 are from Cycle 1 1 proposal ID 9290 and 
CL J1226 + 3332 is from Cycle 10 proposal ID 9033 and are down- 
loaded from HST archive in reduced form. The reason we choose to 
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analyse HST data rather than ground based observations is because 
the affects of atmospheric seeing are removed with only the rela- 
tively small instrumental point spread function (PSF) to consider, 
which makes fitting the profile shape considerably easier at small 
radii and with no near-IR sky background to contend with these 
HST observations are thus ideal for this project. 

We define a low redshift comparison sample of 19 clusters 
at 2median = 0.23 for which archival HST/ACS data are avail- 
able (Table [2}. These are a subset of the clusters select ed for study 



by th e Local Cluster Substructure Survey (LoCuSS, Smith et al 



2010) from t he ROSAT All-sky Su rvey catalogues dEbeling et al 



19981 l200d : iBohringer etal] |2004) . These clusters were all ob- 
served with HST/ACS through the F606W filter which gives good 
rest frame filter agreement with our high redshift F850LP sample. 
The majority of these data were acquired in Cycle 15 (PID:10881; 
P.I.: G. P. Smith) with the balance from Cycle 11 (A 611, Z2701, 
A 2537; PID:9270; P.I.: S. W. Allen) and Cycle 13 (PID: 10200; 
1ES0657-558; PL: C. Jones). All of the observations employed 
standard dither patterns, and spanned at least 1200sec. The flat- 
fielded frames were reduced using standard STScI MultiDrizzle 
routines (Hamilton-Morris et al., in prep.). We choose to anal- 
yse our own low redshift sample rather than compare to those in 
the literature due to the confusion about methods of effective ra- 
dius measurement and the various resulting discrepancies seen be- 
tween studies at low- z discussed in the introduction. The average 
X-ray luminosity of the low redshift clusters in the 0.1 — 2.4 keV 
band is ~ 7 x 10 44 ergs _1 with an average X-ray temperature of 
6.5 ± 2.2 keV which is similar to that of their high redshift counter- 
parts (average L x ~ 6x 10 44 ergs _1 andTx = 7.2±1.8 keV) and 
as these properti es correlate with cluster mass, albeit with a signifi- 
cant scatter (e.g.[R eiprich & Bohringer 2002; Ma ughan et al . 2006; 
IVikhlinin et alj20 09 ; Hoek stra et al.l201 ll) . these can be thought of 
as well matched samples. 



2.1 BCG selection 

The identification of BCGs within clusters is usually obvious from 
visual inspection of the images since, for such rich clusters, they 
are the prominent galaxy closest to the X-ray centroid, often with a 
cD-like profile. However for our high redshift sample we chose to 
formalize this by studying the tip of the red sequ ence in the colour- 
magnitude relations given in IStott et alj feOlCl) . For each cluster 
we identified the red sequence with J — K s colour and selected 
the brightest galaxy from the if s -band magnitudes of all the red 
sequence galaxies wit hin a projected dista nce of 500 kpc from the 
cluster X-ray centroid. lLin & MohJ ( l2004l) have shown that approx- 
imately 95% of massive clusters the BCG lies within this radius. 
Furthermore these are all spectroscopically confirmed members of 
their host clusters. The BCG selection is more straightforward at 
low-redshift, because the ACS observations are centered on the 
spectroscopically confirmed BCG, many of which lie adjacent to 
strongly-lensed background galaxies. 



3 SURFACE BRIGHTNESS PROFILES 
3.1 Sersic profile 

The su rface photom etry of galaxies is often described by a Sersic 
profile dSersidl968h . 



where I(r) is the intensity, r is the radius from the centre of 
the galaxy, r e is the scale radius, I e is the intensity at r e , n in the 
exponent is a free parameter widely known as the Sersic index and 
b n — 2n — 0.327; a coefficient chosen so that r e is the half-light 
radius defined as the radius which encircles half the light from the 
galaxy dGraham etaill 19961) . 

This can also be written in terms of the surface brightness, 
n(r) as 



(2) 



7(r) = I e exp < -b n 



l/n 



(1) 



where is the surface brightness at r e and c n = 
2.5 6„/ln(10). 

Th e Sersic profile i s a generalisation of the n = 4 case first 
used bv lde Vaucouleursl d 19481) to describe the light profiles of el- 
liptical galaxies and the bulge components of disc galaxies. De- 
pending on the value of n, galaxies can be described as disc-like 
with an n = 1 exponential profile or bulge-like with higher n val- 
ues, where ellipticals are expected to have n — 4 de Vaucouleurs 
profiles. BCGs with cD-like morphologies (i.e. those with emis- 
sion from extended halos) are often best fit by a Sersic profile with 
n > 4 or in some cases by 2 profiles, one for th e central bulge and 
one for the outer envelope (Gonzalez e t alj2005l) , to account for the 
excess light above that expected from n = 4 at large radii. 



4 ANALYSIS AND RESULTS 
4.1 Profile fitting 

Due to the degeneracies between n and r e and the potential for a 
dependency on the surface brightness limit of the images we imple- 
ment our own 1-D fitting process for which we have greater control 
over the input data and par ameters compare d to using 'black box' 
software such as GALFIT dPeng et al. 2002). Having said that we 
first perform 2-dimensional Sersic profile fits using the GALFIT 
(version 3) software package in order to obtain the position an- 
gle and ellipticity, of the BCGs. This software requires reasonable 
initial input parameters such as position, apparent magnitude and 
ellipticit y all of which are estimat ed by first running the SExtractor 
package dBertin & Arnoutslll996l) so that the iterative fitting pro- 
cess converges to the correct solution in the shortest possible time. 

To create 1-D profiles for our individual BCGs we use the 
IRAF STScI package Ellipse to extract isophotal intensities with 
radius, including the ellipticity and position angle parameters de- 
rived from GALFIT, while masking out any nearby galaxies inter- 
actively. We then fit the resulting output 1-D surface brightness pro- 
files along the semi-major axis of the BCG with both general Sersic 
and de Vaucouleur profiles using a least squares fitting routine. To 
remove the effect of the PSF we only commence fitting beyond 3a 
of the HST PSF (with the average PSF FWHM for our data ~ 0.11 
arcseconds), however we confirm that the results are robust to be- 
ginning this fit further out (> 5a). These fits are performed to a 
consistent surface brightness limit ^ifssolp = 27.0 at z = 1 which 
corresponds to /ipeoew = 25.6 at z — 0.23 when taking into ac- 
count surface brightness dimming, a small k correction due to the 
slight mismatch between t he 'rest-frame' filte r s and a correction for 
passive evolution using a lBruzual & Charlol d2003l) simple stellar 
population (SSP) model with Zf = 3 and solar metallicity, which 
is a good approximation to the stella r component of these galaxies 
dCollins et al.l2009l ; lstott et alj20l(it) . The resulting individual sur- 
face brightness profiles for both the low and high redshift samples 
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Table 1. The high redshift cluster sample. The radii quoted are along the semi-major axis, b/a is the axis ratio of the BCG. 



Cluster 


R.A. 




Dec. 




z 


r e (n = 4) 


r e (free) 


n 


b/a 






(J2000) 






(kpc) 


(kpc) 






CL J0152.7 - 1357 


01 52 


41.0 


-13 


57 45 


0.83 


27.2+0.4 


16.8±0.5 


2.9±0.1 


0.90 


RDCS J0910 + 5422 


09 10 


44.9 


54 22 09 


1.11 


14.4±0.6 


8.5±0.5 


2.5+0.2 


0.73 


MS 1054.4 - 0321 


10 57 


00.2 


-03 


37 27 


0.82 


24.2±0.3 


62.6±5.2 


6.7±0.2 


0.63 


CL J1226 + 3332 


12 26 


58.0 


+33 


32 54 


0.89 


52.6±0.5 


59.4±2.0 


4.3±0.1 


0.70 


RDCS J1252.9 - 2927 


12 52 


54.4 


-29 


27 17 


1.24 


29.7±0.5 


80.3±11.9 


6.8±0.4 


0.80 



Table 2. The low redshift comparison sample. The radii quoted are along the semi-major axis, b/a is the axis ratio of the BCG. 



Cluster 




RA. 


Dec. 
(J2000) 






z 




r e (n = 
(kpc) 


4) 


r e (free) 
(kpc) 


n 


b/a 


Abell2813 




00 


43 


25 


.5 


-20 


37 


02 





292 


87.0+1 


.3 


52.5+1.9 


3.1+0.1 


0. 


59 


Abell 141 




01 


05 


37 


.6 


-24 


40 


49 


0. 


230 


55.3+0. 


.4 


34.6+0.7 


3.2+0.1 


0. 


65 


Abell 2895 




01 


18 


11 


.5 


-26 


58 


11 


0. 


227 


38.8+0. 


.2 


65.8+1.8 


4.9+0.1 


0. 


76 


RXCJ0220.9 


- 3829 


02 


20 


56 


.5 


-38 


28 


54 


0. 


229 


35.8+0. 


3 


30.8+0.6 


3.7+0.1 


0. 


74 


Abell 368 




02 


37 


28 


.2 


-26 


30 


29 


0. 


22 


77.1+0. 


8 


129.4+8.9 


4.9+0.1 


0. 


52 


RXCJ0304.1 


- 3656 


03 


04 


04 


.4 


-36 


56 


28 


0. 


219 


24.6+0. 


1 


95.4+4.3 


7.3+0.1 


0. 


86 


RXCJ0336.3 


- 4037 


03 


36 


16 


.4 


-40 


37 


44 


0. 


173 


12.1+0. 


1 


376.5+72.8 


12.9+0.5 


0. 


95 


RXCJ0449.9 


- 4440 


04 


49 


56 


.9 


-44 


40 


21 


0. 


1501 


57.9+0. 


.3 


40.8+0.8 


3.3+0.1 


0. 


67 


RXCJ0528.2 


- 2942 


05 


28 


15 


.4 


-29 


42 


59 


0. 


158 


60.7+0. 


4 


25.8+0.3 


2.6+0.1 


0. 


59 


Abell 3364 




05 


47 


38 


.1 


-31 


52 


20 


0. 


148 


15.7+0. 


1 


45.4+1.7 


6.4+0.1 


0. 


95 


RXCJ0638.7 


- 5358 


06 


38 


45 


.8 


-53 


58 


23 


0. 


227 


30.1+0. 


.2 


95.9+5.2 


6.9+0.1 


0. 


72 


GC065819 




06 


58 


15 


.9 


-55 


56 


36 


0. 


297 


54.1+0. 


5 


254.6+25.2 


6.6+0.2 


0. 


71 


GC065832 




06 


58 


38 


.1 


-55 


57 


25 


0. 


296 


19.8+0 


.2 


— ¥ oo 


41.0+6.0 


0. 


73 


Abell 611 




08 


00 


56 


.8 


+36 


03 


24 


0. 


288 


66.5+0. 


.3 


27.7+0.3 


2.5+0.1 


0. 


70 


Abell 781 




09 


20 


26 


.8 


+30 


30 


38 


0. 


29 


27.4+0. 


2 


24.3+0.8 


3.8+0.1 


0. 


88 


Zwicky2701 




09 


52 


49 


.1 


+51 


53 


06 


0. 


214 


43.8+0. 


3 


34.4+0.5 


3.6+0.1 


0. 


68 


RXJ 1000.5 4 


-4409 


10 


00 


31 


.3 


+44 


08 


48 


0. 


154 


12.4+0. 


1 


123.8+13.6 


10.0+0.3 


0. 


81 


Abell 2 187 




16 


24 


14 


.0 


+41 


14 


31 


0. 


183 


48.1+0. 


.3 


38.0+0.5 


3.6+0.1 


0. 


68 


Abell 2537 




23 


08 


22 


.3 


-02 


11 


32 


0. 


295 


56.7+0. 


.3 


31.9+0.5 


3.1+0.1 


0. 


73 



are found in the Appendix and the individual r e values for the gen- 
eral Sersic and de Vaucouleur fits are given in Tables [T] and [2] We 
note that these 1-D fits are in agreement with the less carefully mea- 
sured initial output from the GALFIT 2-D fitting software. As with 
previous stu dies we find that Ser sic index and r e are positively cor- 
related (e.g. iGraham et al] [l996). Fitting a surface brightness pro- 
file for a BCG can be very difficult particularly when a dominant cD 
halo component is present as it can sometimes cause the best Sersic 
fit to tend towards n — Y oo, which in practice means n » 10. As 
mentioned in the introduction, the surface brightness profile will 
then resemble a power law of index ~ — 2 as intensity falls with 
radius with the total integrated light diver ging, resulting in very 
large values for r e , as seen by others (e.g. iGonzalez et al]|2005t) 
and here for a number of the low-z clusters. The average semi- 
major axis scale sizes calculated with a biweight estim ator appro- 
priate for non-Gaussian distributions (Beers et al. 1990) are as fol- 
lows: for the high redshift BCGs r e = 26.9 ± 2.3 kpc (n = 4) and 
r e = 57.3 ± 15.7 kpc, n = 4.3 ± 0.9 (Sersic); and for the low red- 
shift sample r e = 43.8+5.4 kpc (n = 4)andr e = 45.5 + 6.9 kpc, 
n = 4.0 ± 0.4 (Sersic). 

To get a robust measurement of the typical size of a galaxy 
at each epoch we also stack our low and high redshift 1-D sur- 
face brightness profiles in order to provide a smoother and deeper 
dataset. We do this first for the high redshift sample, correcting all 
of the surface brightness profiles to z = 1, accounting for the effect 



of cosmic dimmin g, angular scale size and bo th k and evolution cor- 
rections based on a lBruzual & Chariot] ( 2003) SSP model with solar 
metallicity and a formation redshift Zf = 3 as above. The same is 
then done for the low redshift clusters correcting to z = 0.23, the 
average redshift of the sample. We then fit the high and low redshift 
surface brightness stacks, again beyond 3er of the HST PSF, using 
the same 1-D fitting code and find that for the high redshift BCGs: 
r e = 32.1 ± 2.5 kpc when n = 4 and r e = 47.6 ± 13.7 kpc 
with n — 5.4 ± 0.9 for the best fit Sersic model; while for the 
low redshift sample: r e = 43.2 ± 1.0 kpc when n = 4 and for a 
free Sersic fit r e = 57.9 ± 4.5 kpc with n — 4.8 ± 0.2, in good 
agreement with the average of the individual 1-D fits. These fits are 
performed to a deeper surface brightness limit, /ifssolp = 28.0 
at z = 1 corresponding to ^iF606W = 26.6 at z = 0.23 and we 
note that if the fits are instead performed to the original shallower 
limits of the individual fits the results are still in agreement. The re- 
sulting stacked surface brightness profiles and their de Vaucouleur 
and Sersic fits are shown in Figures[T]and[3]with the corresponding 
residuals about these fits shown in Fi gures[2]and[4l We also look at 
the Petrosian radius jPetrosianll 19761) . which has the advantage of 
not being so affected by the background and p hotometric uncertain - 
ties as r e . Fixing 77 = 1.5 for comparison with B rough et alj d2005l) . 
gives r pet = 41.9 kpc and r pe t = 61.7 kpc for the high and low 
redshift stacks respectively (see Figure[5}. Finally to get a further 
non-parametric estimate of the scale size we integrate the total light 
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Figure 1. Surface Brightness profile stack for the high redshift sample with 
the de Vaucouleur and Sersic profile fits plotted (solid and dashed lines 
respectively). 



Figure 2. The difference between the observed surface brightness profile 
and both the de Vaucouleur and Sersic fits for the high redshift stack (circles 
and squares respectively). The vertical dashed line represents 3cr of the HST 
PSF beyond which the fitting commences. 



within the stacks down to the same surface brightness limits as the 
fits and find that the half-light radius ri/ 2 = 23.6 ± 9.6 kpc for the 
high z data and rw2 = 27.0 ± 3.3 kpc for the low-z sample. 

The results of the average scale size, stack fitting and non- 
parametric analysis are found in Table[3] From these results we can 
see that there is evidence for at most a weak trend in scale size 
or profile shape for BCGs over a period of 5 Gyr, between z = 
1 and z — 0.2. If we consider just the de Vaucouleur fits to the 
stacked clusters and the average of the individual de Vaucouleur 
fits, then the high redshift galaxies at z = 1 are still only found 
to be a factor 0.6 — 0.7 the size of their 2 = 0.2 counterparts. 
This is the maximum evolution implied by our data as the Sersic 
fits and integrated light size estimator suggest this is considerably 
less. Clearly there is no evidence for the large growth factors of 
order 4 reported for other early-type galaxies over a similar epochs. 
In addition to this 3 of the individual Sersic fits at high redshift 
show significant cD halos as their n values are significantly greater 
than the n — 4 de Vaucouleur prediction appropriate to a regular 
elliptical galaxy. 



4.2 Kormendy relation 

The Kormendy relation is a slice through the fundamental plane 
that shows how the scale si ze of galaxies depends on the surface 
brightness dKormendvlll977t) . Formally it is a plot of the surface 
brightness at the effective radius (p e ) vs r e for a de Vaucouleur 
surface brightness fit. BCGs are found to lie off the relation for 
normal elliptical ga l axies as they po s sess large extended halo s 
jHoesselet al1ll987l : ISchombed[l987l ; lOegerle & Hoessellll99lh . 
A similar relation to this is BCG size versus luminosity, which for 
BC Gs is found to be steeper than that o f regular elliptical galax- 
ies teernard i et al. 20Q7|; iLiu et al]|2008l) . although, as mentioned 
in the introduction, this is a matter for some debate as authors us- 
ing essentially the same da ta find a significantly shallower slope 
Avon der Linden e t al. 2007). In Figure[6]we plot the Kormendy re- 
lation for our samples. The slope of this relation is parameterised 
as = A + B log r c and for the LoCuSS sample A — 20.3 ± 0.9 
and B — 2.7 ± 0.8, with the errors reflecting the intrinsic scat- 
ter in the relation, not just the formal error of the fit. This is in 
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Figure 3. Surface brightness profile stack for the LoCuSS sample with the 
de Vaucouleur and Sersic profile fits plotted (solid and dashed lines respec- 
tively). 



good agreement with the work of iBrough et alj J2005h who find 
B = 2.60 ± 0.03 but we cannot rule out larger valu es of B found 
in oth er st udies, e.g. B = 3.1 ± 0.1 found by both lHoessel et al] 
Jl987l) and Oeaerle & Hoessell ( fT99ll) and B = 3.44 ± 0.13 found 
bv lBildfell et al.ld2008l) . We include our fit (solid line) and a das hed 
line with the same slope as that found by Bild fell et alj d2008h for 
comparison. It is obviously not possible to perform a fit to the high 
redshift BCGs and look for an evolution in the Kormendy relation 
with only a handful of points but we include them on the plot after 
correc ting for cosmic dimming and k and evolution correction us- 
ing a lBruzual & Charlotl d2003h SSP model with solar metallicity 
and a formation redshift Zf — 3 as above. As expected from the 
surface brightness fitting results, the high redshift points are found 
to occupy the same region of the plot as the low redshift relation. 
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Table 3. Summary of measured parameters. r pet is the Petrosian radius where the Petrosian parameter rj = 1.5. T"i /a i s a non-parametric half-light radius 
derived from the integrated flux. 
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Figure 4. The difference between the observed surface brightness profile 
and both the de Vaucouleur and Sersic fits for the LoCuSS stack (circles 
and squares respectively). 



Figure 6. The Kormendy relation for our sample. The LoCuSS and the 
high redshift BCGs (accounting for cosmological dimming and k and evo- 
lution correction to match the LoCuSS sample) are represented by circles 
and squares respectively. 



2.0 




Figure 5. Petrosian r\ dependency with radius for the low and high redshift 
stacks (circle and square points respectively). The horizontal line is for r\ = 
1.5 where we measure the Petrosian radius. 



5 DISCUSSION 

The results of our investigation demonstrate that when the best 
available imaging data are used there is evidence for at most only 
a small increase in the scale size of BCGs over the 5 Gyrs between 
z — 1 and z = 0.2 and no evidence for a change in the shape of 
their light profiles. This is in contrast with other BCG studies that 



probe the scale siz e to moderate redshifts (e.g. lNelson et alj|2002l ; 
lAscaso et al. 2010) and with the res ults for massive field galaxies 
(e.g. van Dokkum et al . 2009, 2010). The interpretation of this re- 
sult when taken in concert with the findings of our near-infrared 
observations, which demonstrate that there has been little evolu- 
tion in the stellar mass over the same period, suggests that major 
merging is unlik ely to be an important process fo r BCG evolution 
since z = 1 dCollins et al"]|2009l: Istott et al.ll2010l) . This contrasts 
with semi-analytic models of BCG evolution based on hierarchical 
N-body simulations which require ~ 70% of t he final BCG stel- 
lar m ass to be accreted over this time frame jPe Lucia & Blaizol 
1 20071) . Other models have tried to predict the scale size evolution 
for massive elliptical gal axies and also find significant growth fac - 
tors since redshift 1 (e.g. lKhochfar & Silkll2006l : lNaab et ai]|2009h . 

As discussed in the introduction, a detailed comparison with 
other results is difficult due to the heterogeneity of the samples and 
analyses, so here we restrict ourselves to comparing general trends . 
Our results do not support the findings of Valenti nuzzi et alj d2010h 
who report an average r e value of only 8 kpc using 8 clusters from 
the ESO Distant Clusters Survey at 0.4 < z < 0.8, however this 
sample is not X-ray selected and may contain BCGs which do not 
trace the more massive cD system s associated with X-ray samples 
at high redshift jBurke et alj2000h . 

There are two other recent res ults reporting e volution of BCG 
size at intermediate r edshift: first. lBernard"Hl2009l find that BCGs 
in SDSS optically selected clusters at z ~ 0.25 are smaller than 
their lower redshif t counterparts by as much as 70%; secondly, 
lAscaso et al] d2010h find a reduction of 0.6 in relative scale size be- 
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tween 20 BCGs in the range 0.3 < z < 0.6) and an X-ray matched 
sample at z ~ 0.05. This would be consistent with our result if the 
ev olution takes pla ce relatively recently (z < 0.25), as suggested 
by Bernardi d2009f) . Having said that we do find on average larger 
BCG r e values at z = 0.2 than in these studies and the average 
r e (44 ± 5 kpc) for our LoCuSS sample is more consistent with 
the average r e value (52 i 4 kpc) found bv lBildfellet"aIU2008l) 
for 48 X-ray luminous clusters at an average redshift of 0.26 and is 
arguably a better sample to compare with LoCuSS. 

Restricting our analysis to only the highest quality data avail- 
able for the highest redshift X-ray clusters, we have found that 
no great size evolution exists for BCGs. As far as the result for 
the field galaxies goes we can only speculate that either there is a 
strong environmental dependence, with massive galaxies in clus- 
ters already being morphologically mature at high redshift whereas 
massive field galaxies are yet to undergo a transformation. On the 
other hand, as we emphasised in the introduction, it is not possible 
to rule out observational bias. 

We note that the extended cD halo is often found to modify 
the BCG profile to have a higher Sersic index than the standard 
elliptical-like de Vaucouleur n — 4. As 3 of the individual Sersic 
fits at high redshift have n > 4, we conclude that the cD halo is 
in place for at least some BCGs at high redshift, although these 
are washed out to some extent in the stack. However, some mod- 
els predict that this halo forms late and results in BCGs only re- 
cently departing from the Korm endy relation for normal ellipticals 
dRuszkowski & Springell 120090 . It would be interesting to obtain 
deep data in other optical bands to see whether this cD halo has a 
significant colour gradient at high redshift compared to local sam- 
ples as then we may begin to investigate its age and origin if signif- 
icantly different in colour to the main bulge of the BCG. 

Due to our clusters being high redshift X-ray emitters discov- 
ered in flux-limited X-ray surveys, there is an argument that they 
represent the most massive and therefore rarest objects at that epoch 
and therefore we may not be making a fair comparison. One may 
expect that these massive clusters would themselves contain the 
most massive and therefore by inference the largest sized BCGs at 
that epoch. However, when we look at our low redshift sample we 
find no correlation be tween X-ray lumin osity or temperature and 
r e in agreement with lAscaso et alj|2010] and so we conclude this 
is not a significant issue. One other aspect is the affect that clus- 
ter cool-core strength has on the properties of BCGs, particularly 
as the number of cool-core clusters appears to evolve ( Sant os et al.l 
2010), possibly resulting in inconsistent X-ray selection with red- 
shift. However it is reassuring that a quick check using 18 of the 
low-z clusters observed with Chandra, shows no correlation be- 
tween the inner gas density profile slope and the BCG scale size 
(A. Sanderson, private communication). Another possible bias is 
that to be discovered by X-ray surveys at such high redshifts, we 
are only seeing the most relaxed clusters with the highest central 
gas densities which may be expected to host more morphologically 
mature BCGs at that epoch. Of course we cannot account for this 
and other biases with such a small sample selected from an un- 
known mass function at high redshift. A larger sample of BCGs 
is required above z = 1 , including those hosted by less relaxed 
cluster systems, to account for potential biases. 

Our evidence supports assertion that BCGs do not change in 
appearance over the last 6Gyr and occupy a similar part of the 
Kormendy relation over this entire period. Interestingly the other 
galaxy type for which the Kormendy relation appears not to have 
changed over a similar timeframe is radio galaxies. Although these 
objects have scale sizes less than BCGs (r e ~ 12 — 15 kpc at 



z ~ 0.5 - 1 lMcLureetaf]|2004l) . only modest size evolution 
has been seen in powerful radio galaxies, gr owing from ~ 8kp c 
at 2 = 2 a factor of 1.5 by the present day dTargett et al.ll2010]) . 
This may be suggestive that radio galaxies and BCGs evolve along 
similar evolutionary paths, at least at late times. 

We have a picture of large BCGs residing in the cores of mas- 
sive clusters with their cD halos already in place by z = 1 and any 
subsequent merger or interaction events having, on average, very 
little effect on the global properties of the galaxy. This may be un- 
derstood if we assume that major merger events in the cores of mas- 
sive clusters atz = 1 will be just as rare as their local counterparts 
due to the high velocity dispersions involved. Significant amounts 
of star formation will also be suppressed by the hot cluster environ- 
ment so we can perhaps think of these systems as already fully ma- 
ture at this epoch. The discovery of galaxies with large stellar mass 
and scale size at these epochs is a major challenge to the current 
theoretical models, although comparisons using deeper imaging of 
larger cluster mass-selected samples are required to make further 
definitive tests. 



ACKNOWLEDGMENTS 

Firstly, we would like to thank the anonymous referee for their use- 
ful comments and suggestions. 

We thank Al Sanderson for providing the cool-core strength 
information mentioned in the discussion. 

We acknowledge financial support from Liverpool John 
Moores University and the STFC. GPS acknowledges support from 
the Royal Society 

IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



REFERENCES 

Ascaso B., Aguerri J. A. L., Varela J., Cava A., Bettoni D., Moles 

M., D'Onofrio M., 2010, ArXiv e-prints 
Beers T. C, Flynn K., Gebhardt K, 1990, Astron. J., 100, 32 
Bernardi M., 2009, MNRAS, 395, 1491 

Bernardi M., Hyde J. B., Sheth R. K, Miller C. J., Nichol R. C, 

2007, Astron. J., 133, 1741 
Bertin E., Amouts S., 1996, A&AS, 117, 393 
Bildfell C, Hoekstra H., Babul A., Mahdavi A., 2008, MNRAS, 

389, 1637 

Bohringer H., Schuecker R, Guzzo L., Collins C. A., Voges W., 
Cruddace R. G., Ortiz-Gil A., Chincarini G, De Grandi S., Edge 
A. C, MacGillivray H. T, Neumann D. M., Schindler S., Shaver 
P., 2004, A&A, 425, 367 

Brough S., Collins C. A., Burke D. J., Lynam P. D., Mann R. G., 
2005, MNRAS, 364, 1354 

Bruzual G., Chariot S., 2003, MNRAS, 344, 1000 

Burke D. J., Collins C. A., Mann R. G., 2000, Astrophys. J. Let., 
532, L105 

Carter D., 1977, MNRAS, 178, 137 

Cimatti A., Cassata P., Pozzetti L., Kurk J., Mignoli M., Renzini 
A., Daddi E., Bolzonella M., Brusa M., Rodighiero G., Dickin- 
son M., Franceschini A., Zamorani G, Berta S., Rosati P., Hall- 
iday C, 2008, A&A, 482, 21 

Collins C. A., Mann R. G, 1998, MNRAS, 297, 128 



8 J. P. Stott et al. 



Collins C. A., Stott J. P., Hilton M, Kay S. T„ Stanford S. A., 
Davidson M., Hosmer M., Hoyle B., Liddle A., Lloyd-Davies 
E., Mann R. G., Mehrtens N., Miller C. J., Nichol R. C, Romer 
A. K., Sahlen M., Viana P. T. P., West M. J., 2009, Nature, 458, 
603 

Daddi E., Renzini A., Pirzkal N, Cimatti A., Malhotra S., Sti- 
avelli M., Xu C, Pasquali A., Rhoads J. E., Brusa M., di Serego 
Alighieri S., Ferguson H. C, Koekemoer A. M., Moustakas 
L. A., Panagia N., Windhorst R. A., 2005, Astrophys. J., 626, 
680 

Damjanov I., McCarthy P. J., Abraham R. C, Glazebrook K., 
Yan H., Mentuch E., Le Borgne D., Savaglio S., Crampton D., 
Murowinski R., Juneau S., Carlberg R. G., J0rgensen I., Roth K., 
Chen H., Marzke R. O., 2009, Astrophys. J., 695, 101 

De Lucia G., Blaizot J., 2007, MNRAS, 375, 2 

de Vaucouleurs G., 1948, Annales d'Astrophysique, 11, 247 

Driver S. P., 2010, ArXiv e-prints 

Ebeling PL, Edge A. C, Allen S. W., Crawford C. S., Fabian A. C, 

Huchra J. P., 2000, MNRAS, 318, 333 
Ebeling H., Edge A. C, Bohringer H., Allen S. W., Crawford 

C. S., Fabian A. C, Voges W., Huchra J. P., 1998, MNRAS, 301, 

881 

Fasano G., Bettoni D., Ascaso B., Tormen G., Poggianti B. M., 
Valentinuzzi T., D'Onofrio M., Fritz J., Moretti A., Omizzolo 
A., Cava A., Moles M., Dressier A., Couch W. J., Kjasrgaard P., 
Varela J., 2010, MNRAS, 404, 1490 

Franx M., van Dokkum P. G., Schreiber N. M. E, Wuyts S., Labbe 
I., Toft S., 2008, Astrophys. J., 688, 770 

Gonzalez A. H., Zabludoff A. I., Zaritsky D., 2005, Astrophys. 
J., 618, 195 

Gonzalez A. H., Zaritsky D., Zabludoff A. I., 2007, Astrophys. 
J., 666, 147 

Graham A., Lauer T. R., Colless M., Postman M., 1996, Astro- 
phys. J., 465, 534 
Hoekstra H., Donahue M., Conselice C. J., McNamara B. R., Voit 

G. M., 201 1, Astrophys. J., 726, 48 

Hoessel J. G., Oegerle W. R., Schneider D. P., 1987, Astron. J., 
94, 1111 

Hopkins P. E, Bundy K., Hernquist L., Wuyts S., Cox T. J., 2010, 
MNRAS, 401, 1099 

Khochfar S., Silk J., 2006, Astrophys. J. Let., 648, L21 

Kormendy J., 1977, Astrophys. J., 218, 333 

Lauer T. R., Faber S. M., Richstone D., Gebhardt K., Tremaine S., 
Postman M., Dressier A., Aller M. C, Filippenko A. V., Green 
R., Ho L. C, Kormendy J., Magorrian J., Pinkney J., 2007, As- 
trophys. J., 662, 808 

Lin Y.-T., Mohr J. J., 2004, Astrophys. J., 617, 879 

Liu F. S., Xia X. Y, Mao S., Wu H., Deng Z. G., 2008, MNRAS, 
385, 23 

Lopez-Cruz O., Yee H. K. C, Brown J. P., Jones C, Forman W., 

1997, Astrophys. J. Let., 475, L97+ 
Mancini C, Daddi E., Renzini A., Salmi E, McCracken H. J., 

Cimatti A., Onodera M., Salvato M., Koekemoer A. M., Aussel 

H. , Le Floc'h E., Willott C, Capak P., 2010, MNRAS, 401, 933 
Maughan B. J., Jones L. R., Ebeling H., Scharf C, 2006, MNRAS, 

365, 509 

McLure R. J., Willott C. J., Jarvis M. J., Rawlings S., Hill G. J., 
Mitchell E., Dunlop J. S., Wold M., 2004, MNRAS, 351, 347 

Naab T., Johansson P. H., Ostriker J. P., 2009, Astrophys. J. Let., 
699, L178 

Nelson A. E., Simard L., Zaritsky D., Dalcanton J. J., Gonzalez 
A. H., 2002, Astrophys. J., 567, 144 



Oegerle W. R., Hoessel J. G., 1991, Astrophys. J., 375, 15 
Oemler Jr. A., 1976, Astrophys. J., 209, 693 
Peng C. Y, Ho L. C, Impey C. D., Rix H., 2002, Astron. J., 124, 
266 

Petrosian V., 1976, Astrophys. J. Let., 209, LI 

Reiprich T. H., Bohringer H., 2002, Astrophys. J., 567, 716 

Richstone D. O., 1976, Astrophys. J., 204, 642 

Ruszkowski M., Springel V., 2009, Astrophys. J., 696, 1094 

Sandage A., 1972, Astrophys. J., 173, 485 

Sandage A., 1976, Astrophys. J., 205, 6 

Santos J. S., Tozzi P., Rosati P., Bohringer H., 2010, A&A, 521, 
A64+ 

Saracco P., Longhetti M., Gargiulo A., 2010, MNRAS, pp LI 15+ 

Schombert J. M., 1986, Astrophys. J. Suppl., 60, 603 

Schombert J. M., 1987, Astrophys. J. Suppl., 64, 643 

Sersic J. L., 1968, Atlas de galaxias australes 

Smith G. P., Khosroshahi H. G., Dariush A., Sanderson A. J. R., 
Ponman T. J., Stott J. P., Haines C. P., Egami E., Stark D. P., 
2010, MNRAS, 409, 169 

Stott J. P., Collins C. A., Sahlen M., Hilton M., Lloyd-Davies E., 
Capozzi D., Hosmer M., Liddle A. R., Mehrtens N, Miller C. J., 
Romer A. K, Stanford S. A., Viana P. T. P., Davidson M., Hoyle 
B„ Kay S. T., Nichol R. C, 2010, Astrophys. J., 718, 23 

Stott J. P., Edge A. C, Smith G. P., Swinbank A. M., Ebeling H., 
2008, MNRAS, 384, 1502 

Targett T. A., Dunlop J. S., McLure R. J., Best P. N, Cirasuolo 
M., Almaini O., 2010, ArXiv e-prints 

Toft S., van Dokkum P., Franx M., Labbe I., Forster Schreiber 
N. M., Wuyts S., Webb T, Rudnick G., Zirm A., Kriek M., van 
der Werf P., Blakeslee J. P., Ulingworth G., Rix H., Papovich C, 
Moorwood A., 2007, Astrophys. J., 671, 285 

Trujillo I., Conselice C. J., Bundy K., Cooper M. C, Eisenhardt 
P., Ellis R. S., 2007, MNRAS, 382, 109 

Trujillo I., Forster Schreiber N. M., Rudnick G., Barden M., Franx 
M., Rix H., Caldwell J. A. R., Mcintosh D. H., et al., 2006, As- 
trophys. J., 650, 18 

Valentinuzzi T, Fritz J., Poggianti B. M., Cava A., Bettoni D., 
Fasano G., D'Onofrio M., Couch W. J., Dressier A., Moles M., 
Moretti A., Omizzolo A., Kjsergaard P., Vanzella E., Varela J., 
2010, Astrophys. J., 712, 226 

Valentinuzzi T, Poggianti B. M., Saglia R. P., Aragon-Salamanca 
A., Simard L., Sanchez-Blazquez P., D'onofrio M., Cava A., 
Couch W. J., Fritz J., Moretti A., Vulcani B., 2010, Astrophys. J. 
Let.,721,L19 

van der Wei A., Holden B. P., Zirm A. W., Franx M., Rettura A., 

Ulingworth G. D., Ford H. C, 2008, Astrophys. J., 688, 48 
van Dokkum P. G., Franx M., Kriek M., Holden B., Ulingworth 

G. D., Magee D., Bouwens R., Marchesini D., Quadri R., Rud- 
nick G., Taylor E. N, Toft S., 2008, Astrophys. J. Let., 677, L5 

van Dokkum P. G., Kriek M., Franx M., 2009, Nature, 460, 717 
van Dokkum P. G., Whitaker K. E., Brammer G., Franx M., Kriek 
M., Labbe I., Marchesini D., Quadri R., Bezanson R., Uling- 
worth G. D., Muzzin A., Rudnick G., Tal T, Wake D., 2010, 
Astrophys. J., 709, 1018 
Vikhlinin A., Burenin R. A., Ebeling H., Forman W. R., Hornstrup 
A., Jones C, Kravtsov A. V, Murray S. S., Nagai D., Quintana 

H. , Voevodkin A., 2009, Astrophys. J., 692, 1033 

von der Linden A., Best P. N, Kauffmann G., White S. D. M., 
2007, MNRAS, 379, 867 

Whiley I. M., Aragon-Salamanca A., De Lucia G., von der Linden 
A., Bamford S. P., Best P., Bremer M. N, Jablonka P., Johnson 
O., Milvang-Jensen B., Noll S., Poggianti B. M., Rudnick G., 



The sizes of Brightest Cluster Galaxies 

Saglia R., White S., Zaritsky D., 2008, MNRAS, 387, 1253 
APPENDIX A: SURFACE BRIGHTNESS PROFILES 



10 J.P. Stott et al. 




Figure Al. The surface brightness profiles for the individual LoCuSS BCGs 1 . de Vaucouleur and Sersic profile fits plotted (solid and dashed lines respectively). 
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Figure A2. The surface brightness profiles for the individual LoCuSS BCGs 2. de Vaucouleur and Sersic profile fits plotted (solid and dashed lines respectively). 
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Figure A3. The surface brightness profiles for the individual LoCuSS BCGs 3. de Vaucouleur and Sersic profile fits plotted (solid and dashed lines respectively). 
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Figure A5. The surface brightness profiles for the individual high redshift BCGs with the de Vaucouleur and Sersic profile fits plotted (solid and dashed lines 
respectively). 



